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Abstract 

Recent experiments indicate that the tendency toward the formation of unidirectional charge density waves ( "stripes" ) is common to 
various underdoped cuprates. We discuss momentum-resolved spectral properties of valence-bond stripes, comparing the situations 
of ideal and short-range stripe order, the latter being relevant for weak and/or disorder-pinned stripes. We find clear signatures of 
ordered stripes, although matrix element effects suppress most shadow band features. With decreasing stripe correlation length, 
stripe signatures are quickly washed out, the only remaining effect being a broadening of antinodal quasiparticles. This insensitivity of 
photoemission to short-range stripe order may be employed to distinguish it from nematic order, e.g. in underdoped YBa2Cu306+5. 
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1. Introduction 

States with broken lattice symmetries are of vital in- 
terest in the field of high-temperature superconducting 
cuprates, for at least two reasons: Such states reflect in- 
teresting physics inherent to doped Mott insulators, and 
they may compete with d-wave superconductivity. One 
established example is the variety of unidirectional charge- 
density-wave (CDW) states, commonly dubbed stripes 
1 1-5], which break both translation and rotation symme- 
tries of the lattice. Another candidate is a nematic phase, 
where only rotation symmetry is broken. The melting of 
a stripe phase can proceed in a single transition or via an 
intermediate nematic phase [1,2]. 

Experimentally, significant progress has been made over 
the last few years in determining the properties of stripe 
states. Originally, static magnetic stripe order was detected 
in neutron scattering in La2-xSr x Cu04 (214) compounds 
[5-7]. Later, the simultaneous presence of static charge or- 
der was proven by inelastic x-ray scattering [8-10]. Stripe 
order was found to be particularly stable near 1/8 doping, 
with a spatial period of four lattice spacings in the charge 
sector. In other cuprate families, scanning tunneling mi- 
croscopy (STM) experiments [11-14] detected spatial mod- 
ulations (on the surface), with a period close to four. This 
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static order is short-ranged, locally breaks the lattice rota- 
tion symmetry, and has the strongest signal on the bonds 
(instead of sites) of the square lattice of Cu atoms [14]. 

On the theory side, charge segregation in doped Mott in- 
sulators had been in fact anticipated [15], and a plethora 
of papers on the problem appeared after the initial ex- 
perimental observations (see Refs. [2-4] for reviews). Re- 
cent experimental results allowed to constrain theoreti- 
cal stripe models. For instance, the magnetic excitation 
spectrum of stripe-ordered Lai 5 / 8 Bai/ 8 Cu04 [19] could be 
nicely described by a simple model of coupled spin ladders, 
pointing toward bond-centered stripes dominated by lo- 
cal singlet formation [16]. The real-space structure of such 
valence-bond stripes has been argued [17] to be consis- 
tent with the STM data on underdoped I^S^CaC^Os+a 
and Ca2- a; Na a ;Cu02Cl2 [14]. A plausible scenario for these 
materials is that stripes would be slowly fluctuating in 
the absence of disorder, but get pinned (and hence static) 
due to impurity effects [2,18]. Importantly, impurities act 
as random field on charge stripes, which leads to glassy 
short-range order. Even in compounds with well-ordered 
stripes, impurities will limit the low-temperature correla- 
tion length and smear out the finite-temperature order- 
ing transition in the charge sector (analogous to the two- 
dimensional random- field Ising model), as indeed observed 
very recently in Eu-doped La2- a; Sr a ;Cu04 [10]. An obvious 
question is whether other compounds like YBa2Cu3 06+<s, 
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where clear-cut stripe signatures have not been observed, 
also display a tendency toward CDW. Theoretically, it has 
been shown that slowly fluctuating (or disordered) charge 
stripes influence the spin excitations and induce a charac- 
teristic hour-glass shaped spectrum [20] . Such spectra have 
been measured in YBa2Cu30g+<5 [21], indicating the pos- 
sibility that undcrdopcd YBa 2 Cu 3 06+5 is close to stripe 
ordering as well. Interestingly, neutron scattering in de- 
twinned YBa2Cu306.45 has detected signatures of spon- 
taneous breaking of lattice rotation symmetry [22] , which 
have a natural interpretation in terms of nematic order. 

Angle-resolved photoemission (ARPES), routinely per- 
formed on various cuprates, has so far not identified unam- 
biguous signatures of stripe order, and early experimental 
data [23] remained controversial [24,25]. In this paper, we 
shall discuss whether ARPES is suited to detect stripe for- 
mation and to distinguish between tendencies toward stripe 
and nematic states. While various papers on this subject 
are in the literature [26-30], the recent advances in com- 
paring theory and experiment put more constraints on phe- 
nomcnological stripe models, see Sec. 2. We shall employ 
a mean-field model of valence-bond stripes, and also in- 
clude the influence of spatial disorder using a proper order- 
parameter field theory. For perfectly ordered stripes, our 
results show multiple bands, gaps, and weak shadow bands, 
in agreement with earlier work, and we identify the most 
prominent features. However, most traces of stripe order 
are difficult to discern if both horizontal and vertical stripes 
contribute to the signal. Spatial stripe fluctuations further 
smear out the signal, such that damping of antinodal quasi- 
particles (QP) may the only visible effect of short-range 
stripes. Even for unidirectional stripes (i.e. in the presence 
of a lattice anisotropy), the breaking of rotation symmetry 
is "weak" . We propose to use this property for an experi- 
mental distinction between symmetry breaking driven by 
stripes and driven by nematic order, applicable e.g. to un- 
derdoped YBa2Cu306+<5- 

2. Experimental input 

Our phcnomcnological approach is guided by experimen- 
tal input. Let us list the most important results and to- 
gether with the assumptions entering our modelling. 

(i) We shall assume that the order has stripe instead 
of checkerboard character, i.e., locally either horizontal or 
vertical stripes occur - this is consistent with both neu- 
tron and STM data, (ii) Over a significant doping range, 
the charge order has a period of approximately four lat- 
tice spacings [8,10,14] , and we shall restrict the calculations 
to ordering wavevectors K x = (±7r/2, 0), K y = (0, ±7r/2). 
(iii) The typical modulation amplitude in the charge sector 
is ±20 . . . 30%, according to both STM data on underdoped 
Bi 2 Sr2CaCu208+,5 and Ca 2 _ a; Na :r Cu02Cl2 [14] (if the con- 
trast in the tunneling asymmetry is interpreted as den- 
sity modulation) as well as x-ray data on La2- 2: Ba 2: Cu04 
[8] . (iv) We assume bond-centered instead of site-centered 



stripes, based on the STM data [14] and on the modelling of 
neutron scattering data [19] on La 2 _ a; Ba a ;Cu04 [16]. (v) We 
assume that the dominant modulations are on the Cu-Cu 
bonds [14]; in a one-band model this translates into modu- 
lations of bond kinetic, magnetic, and pairing energies (all 
being invariant under spin rotation and time reversal). In 
such a case, the bond modulations will have a d-wavc-likc 
form factor [17,31]. (vi) We assume that dimerization and 
bond order are the driving forces behind stripe ordering 
[32], whereas magnetic long-range order is less important. 
For simplicity, we shall model stripe states without static 
magnetism. Note that this does not mean that we ignore 
local-moment physics, but instead we assume that those 
moments form singlet valence bonds, which we account for 
by modulated hopping within our mean-field theory. Addi- 
tional spin order will only have a weak influence on our re- 
sults (on the mean-field level), as estimated in Sec. 3 below. 

Finally, in-plane anisotropics need to be discussed. In 
compounds without anisotropy between the a and b axes 
of the Cu02 planes, horizontal and vertical stripes are 
energetically equivalent. Then, domain formation is ex- 
pected, and ARPES will average over both stripe direc- 
tions. Exceptions occur in detwinned YBa2Cu306+a (due 
to the presence of CuO chains) and in the LTT phase of 
214-compounds, i.e., in La2- K Ba 2; Cu04and Nd- and Eu- 
doped La 2 - a; Sr a ;Cu04. In the latter cases, stable stripe or- 
der occurs, with correlation lengths of order 100 unit cells 
[10]. However, the preferred stripe direction alternates from 
layer to layer, and an experiment probing multiple layers 
will effectively again average over both stripe directions. 
In detwinned YBa2Cu306+a, the symmetry is orthorhom- 
bic, which offers the possibility to study bulk properties of 
unidirectionally ordered phases. 

3. Mean-field stripe model 

We employ a mean-field model for striped superconduc- 
tors, which is combined with an order-parameter field the- 
ory for the collective CDW degrees of freedom. Both model 
and methodology have been discussed in Ref. [17], and will 
be sketched here only briefly. 

Single-particle properties are calculated from a BCS 
model of fcrmions on the square lattice of Cu atoms: 

S c = j dT J2 [ 5 kA d r+ e k-^) C ka + A k( C k1 C -kl + C ' C -)_ 

k 

(1) 

where summation over spin indices a is implied. The single- 
particle dispersion consists of hopping to first (i), second 
(f), and third (t") neighbors, and fi is the chemical po- 
tential. The pairing is of d-wave type, = Ao^ cos k x + 
Ao y cosfc y with Aoa; = — Aqj, = Aq. 

CDW order is represented by two order-parameter fields 
^x,y{P, T ) f° r horizontal and vertical stripes at wavevectors 
K x , y , such that the real field Q x (r) — Reip x (r)e tKx ' r (sim- 
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Fig. 1. Schematic real-space structure of a valence-bond stripe state 
with a 4 X 1 unit cell. Cu lattice sites are shown as circles, with 
their size representing the on-site hole densities. The line strengths 
indicate the amplitude of bond variables like kinetic and magnetic 
energies. Our mean-field model, Eq. (2), implements spatial on-site 
variations of the chemical potential (via Ki) and spatial variations 
of nearest-neighbor hopping and pairing amplitudes (via K2...5). 



ilarly for Q y ) measures the modulation of both the charge 
density and bond order (i.e., kinetic energy or pairing am- 
plitude), for r on sites and bonds, respectively. The cou- 
pling between fermions and the collective CDW fields Q x>y 
reads 
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The coupling constants K1...5 decide about the electronic 
struture of the CDW state. As stated above, we assume 
bond-dominated stripe order with a d-wave-like form fac- 
tor. Those are induced by K2...5, with the <i-wave character 
encoded e.g. in k 2 = —K3. 

Perfectly ordered static stripes correspond to ip x = 
const, ip y — or vice versa. The commensurate period-4 
situation translates into an 8 x 8 matrix Hamiltonian (4x4 
in the absence of pairing). Before continuing, let us esti- 
mate modulation amplitudes. To obtain a kinetic-energy 
modulation of ±25%, a 15 . . . 20% modulation of the hop- 
ping t is required. To compare with on-site modulations, 
this number needs to be multiplied by the number of neigh- 
bors, z = 4. For t = 0.15 eV, this results in zSt « 0.1 eV. 
For comparison, static magnetism causes a spin-dependent 
on-site (mean-field) potential of ±zJ(S)/2, where J is 
the exchange constant between the spin- 1/2 moments. 
The ordered moment, (S), in the stripe phases of the 214 
compounds has been found by /iSR to be (at maximum) 
half of the moment of the undoped antifcrromagnet [2], 
(S) < 0.15. With J = 0.1 eV, we obtain an upper bound 
for the triplet-channel modulation of 0.03 eV (which would 
be further reduced by stripe disorder and the resulting 
frustration). This is significantly smaller than the modu- 
lation in the singlet channel, providing some justification 
for neglecting static magnetism. 



Fig. 2. Constant-energy cuts at the Fermi level through the spectral 
function A(k,u>) for ideal period-4 valence-bond stripes, with differ- 
ent chemical potentials leading to Fermi surfaces with non-nested 
(left) or nested (right) antinodal pieces, a) Vertical stripes, momen- 
tum resolution 7r/256, energy resolution 2 meV, logarithmic inten- 
sity scale, b) Same as a), but momentum resolution 7r/32 and energy 
resolution 7 meV. c) Same as b), but horizontal and vertical stripes 
are added, d) Same as c), but with a linear intensity scale. A promi- 
nent stripe signature are straight Fermi surface pieces perpendicular 
to the stripe direction near the antinodal points. 

4. Ordered stripes: Electronic spectra 

To set the stage, we start with results for perfectly or- 
dered stripes. We choose band structure parameters t = 
0.15 eV, t' = — 1/4, t" = t/12. Stripes are characterized 
by V>* = (1 + t)/V2, ki = 0, k 2 = -K 3 = 0.03 eV, giv- 
ing a hopping modulation St of ±0.021 eV. Such ideally 
ordered stripes may be expected at temperatures much be- 
low the charge ordering temperature in compounds with 
little quenched disorder and/or strong commensurate lat- 
tice pinning. Candidates are La2-xBa x Cu04 and Eu-doped 
La2- x Sr x Cu04 near 1/8 hole doping. 

Fig. 2 shows Fermi surface (FS) cuts of the spectral func- 
tion A(k,u> = 0) in the normal state (Ao = 0). The left panel 
has /1 = —0.12 eV, resulting in a hole doping of roughly 
0.12, whereas the right panel has /i = —0.02 eV, illustrat- 
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Fig. 3. Momentum-space cuts of A(k,u>) along horizontal (upper panel) and vertical (lower panel) lines in the Brillouin zone for ideal vertical 
stripes, with parameters as in Fig. 2a, left panel. A spectral weight transfer between different subbands is clearly visible in the lower panel, 
i.e. for cuts parallel to the stripe direction. Note the hopping matrix element, t = 0.15 eV, controls the energy scale on the vertical axis; 
larger values may be appropriate to account for the band structure on the scale of 1 eV. 



ing the effect of "nested" antinodal FS pieces (i.e. sepa- 
rated by K x ). Parts a)-d) represent the same data, but 
differently: a) displays the signal for vertical stripes with 
high energy and momentum resolution and a logarithmic 
intensity color scale, b) has reduced resolution, c) shows 
the signal of horizontal and vertical stripes superimposed, 
and d) shows the same on a linear intensity scale. This com- 
parison illustrates that the presence of both horizontal and 
vertical stripes makes the identification of stripe signatures 
difficult, and that shadow band features are present, but 
generically weak due to matrix element effects. 

In Fig. 3, we display A(k,ui) as function of energy along 
horizontal and vertical cuts in momentum space, for verti- 
cal stripes with parameters as in Fig. 2a left. Finally, Fig. 4 
shows constant-energy cuts through A{k, w) at different en- 
ergies below the Fermi level, but now for superconducting 
stripes with Ao = 24 meV, — K4 = K5 = 7.5 meV. 

Figs. 2, 3, and 4 nicely illustrate all ARPES features 
arising from stripes. Those are: (i) broken rotation sym- 
metry in momentum space, (ii) multiple bands and gaps, 
arising from Bragg reflection at the boundaries of the re- 
duced Brillouin zone, with the gaps being smallest along 
the momentum-space diagonal due to the d- wave-like form 
factor, (iii) shadow bands shifted by the ordering wavevec- 
tor, which are weak in most parts of momentum space, (iv) 
straight FS segments appearing near the antinodes in the 
direction perpendicular to the stripes, (v) nodal QP sur- 
vive at the FS. Thus, the presence of nodal QP cannot be 
taken as evidence against stripes. Closed FS pockets, which 
exist for certain parameter combinations, are essentially in- 
visible due to the small intensity. Antinodal FS pieces in 
the direction parallel to the stripes may disappear in situ- 
ations close to nesting, due to a Bragg gap opening. Then, 
adding horizontal and vertical stripe signals leads to arc- 
like FS pieces (Fig. 2 right). Otherwise, the antinodal FS 
pieces display rather little stripe-induced shift in momen- 



tum space: The Fermi-momentum locations along (0,7r)- 
(7r,7r) and (tt, 0)-(ir, 7r), determined from the most intense 
bands in Fig. 2a left, differ by only Akp w 7r/20. Signatures 
of broken rotation symmetry are more pronounced e.g. at 
0.15 eV below the Fermi level, with an overall oval-shaped 
constant-energy contour. 

Most features discussed here are qualitatively similar to 
earlier work [26-29] . The visible differences between bond- 
and site-centered stripes are minimal, being restricted to 
some matrix elements. Accordingly, our FS structure is re- 
lated to that found in recent work [30] geared towards un- 
derstanding of the quantum oscillation experiments, with 
the difference that our treatment does not include static 
magnetism (which is important in forming small FS pock- 
ets, but does no drastically influence the overall ARPES 
signal as discussed above). 
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Fig. 4. Constant-energy cuts at different energies through A(k,w) of 
ideal superconducting stripes with parameters as in Fig. 2 (left), but 
Aq = 24 meV, — K4 = K5 = 7.5 meV. 
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In our subjective judgement, the most prominent stripe 
features, possibly detectable experimentally, are the 
straight FS segments near the antinodes (Figs. 2) and the 
spectral weight transfer between the subbands as function 
of k x (Fig. 3). The latter feature implies that, upon chang- 
ing k x , different high- intensity subbands cross the Fermi 
level - these bands constitute the near-nodal FS pieces 
and the straight antinodal pieces, respectively. 

5. Fluctuating stripes 

Modelling disordered stripes requires to account for fluc- 
tuations oiip x ^y. Those we assume to be described by a tp 4 - 
type theory for the 0(4) field tJj = (ip x ,ip y ) [18,20,34]. 
The precise form of will determine the character of the 
fluctuations (amplitude vs. phase). As argued elsewhere 
[17,20], a reasonable assumption is that amplitude fluctua- 
tions of ip are small, i.e., stripe disordering is dominated by 
dislocations and domain walls. Hence, we resort to a form 
of with suppressed amplitude fluctuations: 



\d T A\ + |3rVvl + + Sy\lpy\ 2 (3) 



S,p = J drd f 
+c\ x \d x ip x \ 2 4 



^ix \ u y 



\dytpx\ 



z iy l^y^y 



'-2}/ 



\d x ip y 



with ip 2 = \ip x \ 2 + \ip y \ 2 - A combination of u\ < and 
U2 > suppresses amplitude fluctuations of ip. The quar- 
tic v|'!/' a ;| 2 |'!/'i/| 2 term regulates the repulsion or attraction 
between horizontal and vertical stripes; we shall employ 
v > leading to stripe-like order (whereas v < results in 
checkerboard structures) . The phase-sensitive w term pro- 
vides commensurate pinning and selects bond-centered (in- 
stead of site-centered) stripes for w > [14,20]. (Small in- 
commensurabilities enhance phase fluctuations, but do not 
qualitatively change our results.) 

6. Disordered stripes: Electronic spectra 

A full treatment of the action S c +S^+S c ^, is computa- 
tionally difficult. As in Refs. [17,20], we adopt the adiabatic 
approximation that the collective CDW field if does not 
fluctuate on the time scale of the c dynamics, being justi- 
fied for disorder-pinned stripes. Practically, lattice Monte 
Carlo (MC) simulations of on lattices with 64 2 sites are 
used to generate configurations of ipx.y, for each configura- 
tion S c +S C ip is diagonalized to obtain the electronic spec- 
trum, which is then averaged over 50 configurations. The 
underlying picture of short-range-ordered static stripes ac- 
counts for the existence of stripe segments, checkerboard 
domain walls etc., but neglects stripe dynamics and inelas- 
tic processes. 

To avoid proliferation of parameters, we do not explicitly 
model the impurity pinning. Instead we simply work in a 
regime of with short-range order, and assume that the 
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Fig. 5. Constant-energy cuts through A(k,u>) of short-range ordered 
non-superconducting stripes, with band parameters as in Fig. 2. The 
differences between stripe correlation lengths of £ ss 50 (left) and 
£ 4 (right) are minimal, and stripe features are barely visible. 
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Fig. 6. Same as in Fig. 5, but for short-range ordered stripes in the 
presence of a uniaxial anisotropy. 
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configuration average is equivalent to a disorder average. 
[We expect that details of the impurity potential only mat- 
ter for short-range (i.e. point-like) impurities [35], whereas 
a smooth impurity potential arising from out-of-plane de- 
fects essentially pins existing domains.] 

Figs. 5 and 6 show ARPES spectra for different £ [33], for 
situations without and with in-plane anisotropy [induced 
by slightly anisotropic masses, s x . y , in Eq. (3)]. Note that 
£ — > oo corresponds to the stripe ordering transition where 
fluctuations are not small. 

Without anisotropy, Fig. 5, stripe signatures are only 
discernible for very large stripe domains. In contrast, for 
correlation lengths of 30 and below, those signatures are 
smeared, and the strong scattering due to the stripe disor- 
der potential leads to rather incoherent antinodal QP. 

The situation is different in the presence of an in-plane 
anisotropy, Fig. 6: Here, stripe signatures (in particular the 
straight FS pieces near the antinodes) are visible, provided 
that the stripes in the dominant direction have a correlation 
length of order 10 or larger. (Both cases in Fig. 6 have 
an amplitude ratio (ipx)/i' l l>y) °f roughly 2:1, i.e., vertical 
stripes clearly dominate.) 

7. Stripes vs. nematics 

In detwinned samples of underdoped YBa2Cu 3 06+(5, sig- 
natures of spontaneous breaking of rotation symmetry have 
been found in neutron scattering [22]. The spin fluctuation 
spectrum has been modelled assuming a d-wave nematic 
instability [36]. However, it is also conceivable that the pri- 
mary cause of rotation symmetry breaking is the tendency 
toward stripe order (which may be preceded by a nematic 
transition). The difference in these concepts is in the driv- 
ing mechanism of the phenomenon. 

We propose that the two situations can be distinguished 
in ARPES. Fig. 2 shows that the horizontal and vertical 
antinodal Fermi wavevectors for "realistic" unidirectional 
stripes differ only by roughly 7r/20. In a purely nematic 
picture, the order parameter is the effective anisotropy in 
the hopping matrix elements t x and t y , distorting the FS. 
The difference between the horizontal and vertical Uf is 
linear in At = (t x — t y )/2, and At/t w 4% is sufficient to 
produce an antinodal Akp of 7r/20. However, in Rcf. [36] 
the effective At/t required to model the neutron data was 
about 20%, which does not only shift the antinodal fc^, but 
even leads to a change in the FS topology. 

Independent of precise numbers, it is clear that ARPES 
will be able to distinguish whether the driving mechanism 
of the anisotropy seen in Ref. [22] is of stripe nature (lead- 
ing to little antinodal shifts) or of purely nematic nature 
(implying large antinodal shifts). 

8. Conclusion 

We have discussed photoemission signatures of ordered 
and statically disordered valence-bond stripes in cuprates, 



using a combination of mean-field theory for the fermionic 
sector and Monte-Carlo simulations for the stripe order pa- 
rameter. The results show clear stripe fingerprints for per- 
fectly ordered stripes, in qualitative agreement with earlier 
work. These fingerprints should be visible in high-resolution 
experiments on stripe-ordered 214 compounds. However, 
moderate stripe disorder, in combination with the simul- 
taneous presence of horizontal and vertical stripe domains, 
is able to wash out most stripe signatures. This implies 
that stripes with a spatial correlation length of e.g. 10 lat- 
tice spacings, which are easy to see in an STM experiment, 
leave little trace in ARPES. We have proposed that the 
position of antinodal quasiparticles can be experimentally 
used to distinguish whether rotation symmetry breaking in 
YBa2Cu306+<5 is driven by stripe or nematic order. 
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